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ABSTRACT 

The unitary NN—nNN model contains a serious theoretical flaw: unitarity is obtained 
at the price of having to use an effective ttNN coupling constant that is smaller than the 
experimental one. This is but one aspect of a more general renormalization problem 
whose origin lies in the truncation of Hilbert space used to derive the equations. Here 
we present a new theoretical approach to the irNN problem where unitary equations are 
obtained without having to truncate Hilbert space. Indeed, the only approximation made 
is the neglect of connected three-body forces. As all possible dressings of one-particle 
propagators and vertices are retained in our model, we overcome the renormalization 
problems inherent in previous 7rNN theories. The key element of our derivation is the 
use of convolution integrals that have enabled us to sum all the possible disconnected 
time-ordered graphs. We also discuss how the convolution method can be extended to 
sum all the time orderings of a connected graph. This has enabled us to calculate the 
fully dressed AW one pion exchange potential. We show how such a calculation can be 
used to estimate the size of the connected three-body forces neglected in the new tiNN 
equations. Early indications are that such forces may be negligible. 

INTRODUCTION 

Eighteen months ago, at the International Few-Body Conference in Adelaide, one of us 
(B.B.) summarized the status of our understanding of the NN — tiNN system as follows 
"Despite the almost thirty years of effort, we are forced to the conclusion that all 
current methods have serious deficiencies and that a radical new approach may be needed 
before this outstanding problem is finally resolved". In particular, it was demonstrated 
that the so called "unitary NN — ttNN model" [2-6], which had been the state of the 
art description for more than 10 years, is theoretically flawed. It is therefore especially 
pleasing to be able to present, to this conference, just such a radically new approach to the 
ttNN system, involving convolution integrals, which resolves the outstanding theoretical 



Figure 1: Allowed dressing in the unitary NN — ttNN model, with associated Z renor- 
malization factors, (a) ttN nucleon pole graph, (b) NN OPE graph. 



problems of the unitary NN—irNN model, and which, at this stage, promises to effectively 
take into account all the diagrams of time-ordered perturbation theory.[| 

As it still may not be well recognized that the unitary NN— kNN model is flawed, let 
us first briefly repeat here the nature of the theoretical error inherent in the model, and 
how this error can lead to large inaccuracies in predictions of observables. 

The origin of the problem lies in the truncation of Hilbert space used to derive the 
NN — ttNN equations. Truncation, in this case, means that explicit states with more 
than one pion are forbidden (certain multi-pion states are, however, included implicitly 
through the use of two-body, energy- independent potentials) . This truncation has serious 
consequences for the renormalization of both the two-nucleon propagator and the ttNN 
vertex. In Fig. 0(a) we show the ttN nucleon pole diagram where the intermediate state 
nucleon is dressed by one-pion loops; however, the initial and final state nucleons do not 
include dressing since two-pion states are neglected in the truncation. Since close to the 
nucleon pole the dressed one-nucleon propagator is of the form g(E) ~ Z/(E — m), where 
Z is the residue at the pole, Fig. 0(a) illustrates how each irNN vertex f(E) gets effectively 
renormalized by a factor of Z 1 ^ 2 . Thus f n NN — Z 1 l 2 f{m) is essentially the -kNN coupling 
constant, and this fact is used to fix the strength parameter in the form factor f(E). 
With all other parameters of f(E) fixed to reproduce experimental irN phase shifts, this 
form factor then enters the unitary NN — ttNN equations as an input. As illustrated in 
Fig. 0(b), when the NN one pion exchange (OPE) amplitude is calculated in the unitary 
NN— TtNN model, the initial and final nucleons are dressed by pions and consequently 
each external nucleon obtains a renormalization factor of Z 1 ^ 2 . The first renormalization 
problem is the fact that Z ^ Z. This arises because two nucleons cannot be dressed at 
the same time in the truncated Hilbert space; thus, each nucleon in a two-nucleon state 
cannot obtain its full dressing. This, however, may not be such a serious problem since, in 
practice, the difference between Z and Z turns out to be quite small. The serious problem, 
instead, is the size of the effective irNN coupling constant in the NN — nNN equations. 
Taking Z m Z, Fig. 0(b) illustrates that each vertex gets renormalized by a factor of Z, so 
that the effective irNN coupling constant here becomes Zf(m); this is a factor Z 1 / 2 times 
the physical coupling constant. With Z being typically between 0.6 and 0.8, we come 
to the disturbing conclusion that the effective irNN coupling constant in the NN— ttNN 
equations is smaller than the one used in constructing the input. This observation helps 
explain why one typically obtains much too small pp — > n + d cross sections using this 
model [8-11]. 

These important observations about the renormalization problem in the unitary NN— 
ttNN model were already made in 1985 by Sauer et al. yet they seem to have 



gone largely unnoticed. Perhaps this is partly because one could find less "fatal" reasons 
for the low cross sections; for example, it was legitimately argued that off-shell effects 
and the lack of a "backward-going pion" in the AW — > iVA amplitude can lead to the 
underestimation of pp — >7i + d cross sections [IU|. However, with the advent of calculations 

1 We note that progress has also been made in the four-dimensional relativistic sector where recently 
derived covariant equations for the nNN system resolve overcounting and undercounting problems of pre- 
vious formulations H. In this review talk, however, we shall limit the discussion to the three-dimensional 
approaches of time-ordered perturbation theory. 



where the nucleoli and A are treated on an equal footing in the NN — irNN model (the 
so-called BB— irBB equations the effective irNN coupling constant is lowered by yet 
a further factor of Z 1 / 2 in the most important NN — >NA amplitude, and it has become 
very apparent that the renormalization problem is indeed fatal to this type of approach 
to the irNN system. 

It may seem that one can fix the renormalization problem "by hand" by strategically 
including extra Z 1 ^ 2 factors in either irNN propagators or irNN form factors, or both. 
But it soon becomes apparent that there is no easy way of doing this without destroying 
the three-body unitarity of the equations. 

Here we report on a completely different formulation of the irNN problem where 
unitary equations are obtained without having to truncate the Hilbert space to some 
maximum number of pions. Consequently, all possible dressings of one-particle propaga- 
tors and vertices are retained in our model. In this way we overcome the renormalization 
problems discussed above. The key element of our derivation is the use of convolution 
integrals that have enabled us to sum all the possible disconnected time-ordered graphs 



T3| . Indeed, we basically make only one approximation in deriving the irNN convolution 
equations, namely, we neglect connected three-body forces. 

As the details of this new approach have already been published |[14|| , here we shall 



give a shortened derivation, retaining only the essential steps. A substantial part of our 
presentation involves an extension of the convolution idea to connected diagrams. In 
particular, we show how any diagram, connected or disconnected, can be calculated so 
that all possible dressing contributions are included. We then use the derived formalism 
to calculate the AW OPE potential where the nucleons are fully dressed. Not only does 
this calculation indicate that it may be important to include dressing into the popular 
one boson exchange models of the NN interaction, but it also provides for us a first check 
on the size of the neglected three-body forces. We find that, in this case, the three-body 
connected force contribution is negligible. 

THE CONVOLUTION nNN EQUATIONS 

We consider a time-ordered perturbation theory of nucleons and pions described by a 
Hamiltonian H. The exact form of H need not be specified. The Green function for the 
irNN^-irNN process is thus defined by 

(p'l.P^Psl^^lPi.Pa.Ps) = (p'i, P2, Pal E+ \ H |Pi, P2, p 3 ) (1) 

where p Q (p^) denote initial (final) momenta; here, as below, a = 1, 2 label the two 
nucleons while a = 3 denotes the pion. Note that we suppress spin-isospin labels in 
order to save on notation, and it is assumed that the nucleons are distinguishable since 
antisymmetrization can be carried out at the end. 

We define the irNN— >■ irNN t-matrix T(E) by the equation 

G(E) = G {E) + G (E)T(E)G (E) (2) 

where Gq(E) is the free irNN propagator with all particles fully dressed (i.e. it is the fully 
disconnected part of G(E)). Note that in this section symbols G and g shall represent 
operators, while in the rest of the paper they are usually numbers; this should, in any 
case, be clear from the context. Eq. @ is thus a relation between operators acting in the 
momentum space of one pion and two nucleons. In a previous work flT3"|, we have shown 



how Gq{E) can be expressed in terms of the individual pion and nucleon propagators by 
a convolution integral. In particular, if we write 

(p a \g a (E)\ Pa ) = (p'J E+ l _ H \ Va ) (3) 



Figure 2: Diagrammatic representation of potentials (a) V±, (b) V3, and (c) V4. The 
open circles represent all possible graphs excluding those that lead to kNN intermediate 
states. Potentials V 2 and V 5 are obtained by interchanging the two nucleons. 



for the dressed nucleon and pion propagators, then Go(E) is given by 

G (E) = (-ttV) / dz 1 dz 2 g l (E-z 1 )g 2 (z l -z 2 )g 3 (z 2 ). (4) 

Here it is understood that each single-particle propagator is an operator acting in its own 
Hilbert space. 

With T(E) being defined as in Eq. (|2|), we can write 

T(E) = V(E) + V(E)G (E)T(E) (5) 

where the potential V(E) is the sum of all possible ttNN- irreducible graphs excluding 
those consisting of fully disconnected ttNN states. At this point we shall neglect the 
connected diagrams of potential V(E); these correspond to connected three-body forces 
and will be considered in the next section. It should be stressed that the connected 
diagrams of V are the only diagrams that are neglected in this model. 

For the following discussion we follow the labelling convention where amplitudes that 
involve the pion interacting with nucleon i are labelled by subscript i, while the AW 
interaction with a pion spectator is labelled by subscript 3. Further, the indices i, j, 
are reserved for nucleons 1 and 2 while the Greek indices a, (3, 7 go from 1 to 3 unless 
otherwise indicated. 

Of special relevance to the following discussion is the work of Stelbovics and Stingl 
(SS) [pl|, who investigated the disconnectedness problem in the ttNN system, and then 
applied it to a model where Hilbert space is chopped at the two-pion level. As discussed 
by SS, all disconnected 3^3 diagrams belong to one of five classes of disconnectedness, 
denoted by S a , characterized by an appropriate momentum-space 5-function. Three of 
these classes, S\, 5 2 , S3, have two of the particles interacting, the third being spectator. 
The class 64 has the pion in initial state being absorbed by nucleon 2 with the final state 
pion being emitted by nucleon 1. Class £5 is the same as $4 but with nucleons 1 and 2 
interchanged. 

According to this classification, the potential V(E) can be written as the sum 

V(E) = £ V a (E) (6) 

a=l 

where V a (E) consists of diagrams of disconnectedness S a . The potentials V a (E) are rep- 
resented diagrammatically in Fig. [| It is worth mentioning that one of the contributions 
to the potential of Fig. 0(c) is the so called Jennings term which has been shown to be 



important for calculations of the tensor polarization £20 i n scattering ||16|| . 

As it stands, Eq. (|) cannot be used directly for calculations as its kernel is discon- 
nected. To derive equations with a compact kernel we proceed by analogy with the case of 
Faddeev equations and eliminate the potentials V a (E) in Eq. (|5]) in favour of completely 
summed contributions of disconnectedness S a . Let us therefore denote by w a the Green 



function consisting of the set of all diagrams, reducible and irreducible, belonging to the 
disconnectedness class S a . We can then define the corresponding amplitudes w a by 

w a = G w a G . (7) 

As J2a w a is the disconnected part of T, one can use Eq. (|) to equate the disconnected 
parts of type 8 a , and thus obtain relations between the V a and the w a . These relations 
can then be used to reexpress Eq. (El) as [TH 



T = Y,T a (8) 

Q = l 

5 

T a = W a + Y^ W a K af 3G Tp (9) 
/3=1 

where K a p is a 5 x 5 matrix with elements Kn = K22 = ^33 = ^51 = ^42 = k>u = ^25 = 
K45 = K54 = 0, all other elements being equal to 1. The first iteration of Eq. (|9|) then 
gives us a compact kernel. 

As our equations and the ones of SS share the same disconnectedness structure, it is 
especially interesting to compare the actual models used. In the case of the SS model, 
the potentials V a were specified by graphs of the lowest possible order. It then turned 
out that their w a , needed to specify the kernel of Eq. (§), are not directly related to 
usual subsystem t-matrices, and must therefore be specified by solving very complicated 



integral equations [17]. By contrast, we take the most complicated possible choice for the 
V a , namely the set of all possible contributing diagrams of field theory. The remarkable 
thing is, that it is just this maximally complicated choice that enables us to express the 
needed w a as simple convolutions of dressed subsystem amplitudes (the usual two-body 
t-matrices t a , the ttNN vertex /, and the dressed one-nucleon propagator g) |13| . As 
these convolution integrals are written for Green function quantities, we shall utilize a 
"tilde" notation, as in Eq. fl7|), to label amplitudes with additional initial and final-state 
propagators. Thus for the irN t-matrix we define 

U{E) = g« Ni {E)U{E)g« Ni {E) (10) 

where g n N is the dressed nN propagator. We note further that U is the connected part 
of the full uNi Green function, i.e. 

(Pi,P3|ti(£)|Pi,P3) = (Pi,Pal g+ _ g lPi,P3)c- (11) 

In a similar way we can define £3. With these definitions, w a can be written directly 
in terms of convolutions of t a and the spectator particle propagators g a . Introducing the 
short-hand notation 

c = a®b (12) 

to mean the convolution integral 

1 f°° 

c(E) = :/ dza(E - z)b(z), (13) 

2,7ii J —00 

we have that (i ^ j) 

Wi = ti® g i ; w 3 = t 3 (8) g 3 . (14) 

It should be noted that ti(E) is the full ttN t-matrix and can therefore be written in 
terms of nucleon pole and non-pole parts: 



t i (E) = f i (E)g i (E)f i (E)+t b i (E) 



(15) 



where t\(E) is the non-pole part, f%(E) is the dressed vertex for N — > ttN, and fi(E) is 
the dressed vertex for ttN-^N. These vertices are explicitly given in terms of fi and ff 

fi = ShNifi9i ) fi = 9ifi9nN z ( 16 ) 

where 

&,P3\kE)\Pi) = (pJ.Psl^TfflP*> ( 17 ) 

<pj|/ 4 ( J E0|p i ,p 3 > = (Pii^r^lp-P3>- (is) 

In a similar way, we can also express W4 and w& directly in terms of the dressed irNN 
vertex using 

w 4 = fi®f 2 ; w 5 = f 2 ®f 1 . (19) 

As the NN channel is hidden in the input terms w a , the relation of Eq. (Q) to the 
full NN amplitude is not clear. Neither is it apparent if in the process of iteration, the 
intermediate NN propagators will obtain their full dressing. For these reasons, we would 
like to recast Eq. (||) into a form that explicitly exposes the NN channel. 

For this purpose, we write Eq. @ in an AGS form 

5 

U a /3 = K a pG Q + 2J KayW^GoU^p. (20) 
7=1 

We note, however, that in the ^-matrix, rows 1 and 5 are identical, as are rows 2 and 4, 
columns 1 and 4, and columns 2 and 5. This means that 

U a 4 = U a i ; U a s = U a 2 ] Uifj = ] U^p = U\f}. (21) 

Thus although Eq. ( p0|) is a set of 5 x 5 coupled equations, it can be reduced to a set of 
3x3 equations. 

We recall that the terms Wi (i = 1, 2) are defined as the set of all possible disconnected 
graphs for pion scattering on nucleon i. In order to expose two-nucleon states, we define 
the Green function wf as the set of all two-nucleon reducible graphs belonging to Wi, and 
write for the corresponding amplitudes 

Wi = w° + wf (22) 

where w 4 ° is two-nucleon irreducible. Since we consider all possible contributions, it is 
clear that one can write 

wf = FiGNNFi (23) 

where (jaw is the fully dressed two-nucleon propagator, Fi and F iy are the fully dressed 
ttNN vertices in the two-nucleon sector. Again we can write these in terms of the one- 
nucleon propagator gi and ttNN vertices fi and fi through the convolution expressions 

G NN = gi®g 2 , (24) 
G FiG NN = fi®gj, (25) 
G NN FiG = fi® 9j . (26) 

In the same way we can separate out the two-nucleon reducible contributions from w 4 
and w§, and with the further convention that wf = 0, we can write generally 

w„ = w° a + <• (27) 



Clearly 



w^ = F l G NN F 2 ; w£ = F 2 G NN F 1 . (28) 

With these definitions we can define the amplitudes U®p corresponding to all possible 
two-nucleon irreducible contributions, namely 

5 

Ulp = n aP G Q l + K ay w°G U°p. (29) 

7=1 

From Eq. fl2§), Eq. (£7|) and Eq. it follows that 

U*p = C + E UlyGotfGoUtf. (30) 

7=1 

Using the explicit forms for given in Eq. ( |23"D and Eqs. (^), and making use of Eq. 
(|21~1 ) to eliminate channels 4 and 5, we obtain for the last term in Eq. (|3(i|) 

5 

UayG w^ GqIIjp = (U® 1 G F 1 + f/° 2 Go-F 2 ) Gatat (FiG Uip + F 2 G U 2 p). (31) 

7=1 

Substituting into Eq. ( P0|) leads directly to the equation for nucleon-nucleon scattering 

Tnn(E) = V NN (E) + V NN (E)G NN (E)T NN (E) (32) 

where 

Vnn(F) = J^F.GoU^GoFj (33) 
Tnn(E) = Y.F&QUijGoFj. (34) 

This derivation shows explicitly that the NN propagator, which we did not have as an 
explicit input in Eq. (|20|) , is indeed the fully dressed Gm which is given in terms of the 
convolution integral of Eq. (|2"4|). Although it may be clear that the amplitude Uij is made 
up of fully dressed component amplitudes, it is noteworthy that the connection to the 
NN scattering amplitude is via the vertex functions F{ and F{ which themselves have all 
possible dressings. 

A feature of our formulation is that the input to our derived equations, Eq. @ or 
equivalently Eq. (p0[), consists just of usual two-body subsystem off-shell amplitudes. This 
is an important aspect for practical calculations and was not the case for the formally 
similar equations of SS. On the other hand, due to the terms w® and the kernel of the 
equations nevertheless is not of the standard type where all interactions are of two-body 
type. We shall now show how one can rewrite the equations to have standard pair-like 
interactions in the kernel but at the expense of introducing an extra dimension in the NN 
sector. 

We start by writing an alternative decomposition of w a to that given in Eq. (^). For 
Wi we convolute ti, as given in Eq. ([15|) , with a dressed spectator nucleon; that is, the 
decomposition 

Wi = + wf (35) 

will be defined via the convolutions 

G w[G = (g-xNfgfg-K^i ® 9j (36) 
G WiG = (g n Nt b gnN)i® 9j- (37) 



We can now write w a = + generally by defining w£ = 0, wf = w±, w$ = W5, 
= W3, W® = 0, and w® = 0. Analogous to Eq. fl36] ) we have that 

G0W4G0 = (g nN fg)i ® (,9f9*Nh, (38) 
G ws G = (g 7TN fg)2 ® (gfg*N)i, (39) 

which are just more explicit versions of Eq. fll9"|) . 

Introducing this decomposition of u> a into Eq. (p0|), leads us to define, by analogy to 
Eq. (H), 

5 

U% = n a pG Q l + K ai w°G U^. (40) 
7=1 

Analogous expressions to Eq. (|30|) and Eq. (|31|) are also easily written. Using a short-hand 
notation, illustrated generally by 

(ab)i z = a x {E - z)b x {E - z) ; {ab) 2z = a 2 {z)b 2 {z), (41) 

we are then led to define a new kind of NN potential and t-matrix depending on extra 
variables: 

V NN {z\z;E) = Y,U9-N)i Z 'Ul(g.Nf) jz (42) 

ij 

T NN (z',z;E) = ^2{fg^N)iz'Uij(g n Nf)jz- (43) 

ij 

In this way the integral equation for NN scattering acquires an extra dimension, 

T NN (z', z; E) = V NN (z', z; E) 
1 c°° 

-— / dz" V NN (z', z"; E)g x {E - z")g 2 (z")T NN (z" , z; E). (44) 
2m J -00 

How T NN (z', z; E) can be identified with the physical NN t-matrix, T NN (E), can be seen 
as follows. From Eq. (|25|), it is easy to show that in the case when the two initial-state 
nucleons are on-energy-shell, i.e. pf/2m + ^\/2m + 2m = E, then 

g nNi (p 2 j2m + m)f t (p 2 t /2m + m)\p lP2 ) = G (E)F l (E)\p 1 p 2 ) (45) 

with a similar equation holding for fig-nNi- It can therefore be seen from Eq. ( f43D that 
if the initial and final nucleons are on-energy-shell, and if we set the z variables to the 
energies of the second nucleon, i.e. z' = p' 2 2 /2m + m, z = p 2 /2m + m, then 

(piP 2 l T Mv(^,^;^)|PiP2) = (pip 2 |T MV (E)|p 1 p 2 ) (46) 

where Tnn(E) is given in Eq. (|34]). In this sense Tnn{z' , z\ E) can be considered as an 
off-shell scattering amplitude with additional energy-like variables z and z' . 

Our final result for the description of the ttNN system consists of the coupled equations 
(|40| ) and (0). The equation for U° is just the Faddeev equation for the nNN system with 
no absorption. In momentum space it is a 6-dimensional integral equation and involves 
pair-interactions n-N and N-N. 

The equation for AW scattering, Eq. (PD, in momentum space, is a 4-dimensional 
integral equation. In this sense it is similar to the Bethe-Salpeter equation but with 
particular forms for the nucleon propagator and NN potential. In the case of Eq. (f44|), 
the nucleon propagator contains only positive energy states while the one-pion exchange 
potential, for example, has irNN vertex functions that depend only on one energy variable 



- in the Bethe-Salpeter case they depend on two energy variables. It is interesting to note 
that by introducing a fourth dimension into the NN sector, we take into account the NN 
irreducible diagrams in w 4 and u>5 while at the same time having only pair-like interactions 
in the kernel for U°. 

We conclude by noting that our derivation is model independent. Thus, for example, 
there was no need to specify either how pions actually couple to nucleons, or what model 
is used for dressing the pion. We also note that our expressions for the coupled %NN 
equations, Eq. fllPf) and Eq. (|44|), or alternatively the coupled equations ( ^9|) and ([34"D , 
express the main features of the physics. Thus, for example, we can easily check that 
our formulation gives consistent renormalization, simply by substituting the pole term 
Z/(E + — m) wherever the single-nucleon propagator g(E) appears. The unitarity of our 
equations is also evident from their derivation. For practical calculations, however, either 
of the coupled set of equations can easily be cast into one set of coupled equations for 
the reactions NN — > NN, nd^-nd, and NN^nrd. In the case of Eq. fl4D|) and Eq. ([0]), 
this would result in coupled equations that are of similar form to the well known unitary 
NN — 7rNN equations [||], but with an extra dimension in the NN sector. 



THE DRESSED NN OPE POTENTIAL 

The essential feature of the ttNN equations, presented above, is that they include 
all possible disconnected 7riV7V-irreducible graphs. This is the first time that all such 
disconnected diagrams have been included, and has been achieved only after preliminary 
work showing that convolution integrals can sum all possible time orderings of disconnected 
graphs of time-ordered perturbation theory |TJ| . 



With all the disconnected diagrams included, the focus of attention turns naturally 
to the question of how large the neglected connected 7riV7V-irreducible graphs are. This is 
not an easy question to answer without a systematic analysis of the whole range of such 
connected three-body force contributions. Here we shall take the first steps in such an 
analysis by examining the most basic ingredient of the NN force, namely, the NN one 
pion exchange potential. 

Within the formalism of the ttNN convolution equations, the OPE potential takes the 
form appearing in Eq. (33) with replaced by Gq 1 (the inhomogeneous term of Eq. 



J)). It will be sufficient to examine just one time ordering of the exchanged pion, and 
in this case, the OPE potential of the ttNN convolution equations is given by 

V° PE = F 1 G F 2 . (47) 

The consequence of neglecting connected three-body forces becomes immediately ap- 



parent upon examination of the dressing contributions to this OPE potential. In Eq. (47), 
each of the terms Fi, Gq, and F 2 , is by itself completely dressed. Thus a simple example 
of a dressing contribution that is included in Eq. (|47| ) is given by Fig. 0(a). On the other 
hand, the topologically similar graph of Fig. ^(b) is not included in the formalism of the 
previous section since it is of the form F2V C F\ where V c is the connected irNN — > irNN 
three-body force. 

In order to see the effect of neglecting dressings of the OPE potential like that of Fig. 
[|(b), it would be useful to have an expression for the fully dressed OPE potential. That 
is, we need the OPE graph where all possible dressing graphs involving a ttNN vertex are 
retained (the pion itself, however, will be assumed to have no dressing contributions). As 
will be shown in the sections below, such an expression can indeed be easily derived by 
extending the convolution idea to connected diagrams. Here we just give the result for 



Figure 3: (a) Example of a dressing diagram included in the convolution ttNN equations, 
(b) Example of a dressing diagram involving a connected three-body force - such diagrams 
are not included in the convolution irNN equations. 



the corresponding OPE Green function (in the cm. system) with momentum conserving 
(^-functions removed: 

V° PE (p',p;E) = (-^-) 2 [ dzdz' g(z',p')f(p',p,z',z)g(z,p)^ 

\ iTllJ J Z' — Z — U)f. 

g(E - z', -p')/(-p', -P, E-z',E- z)g(E - z, -p). (48) 

Here u>}~ is the energy of the exchanged pion, and g(E, p) is the dressed nucleon propagator 
defined as the matrix element of Eq. (|3]), but with the momentum conserving (^-function 
removed, i.e. 

S(p-p')g(E,p) = (p'\^-^\p). (49) 

A novel feature of Eq. (^) is the need for a irNN vertex function f(p', p, E', E) that 
depends on two energy variables. By contrast, in the irNN convolution equations, as in 
most other formulations utilizing time-ordered perturbation theory, the dressed vertex 
depends only on one energy variable. The relation between one- and two-energy vertices 
will be given below. 

The essential point here is that if one neglects connected three-body forces from Eq. 
(48|), then one will obtain the OPE Green function of the convolution ttNN equations, 



which numerically is given as a product of three convolution integrals: 

1 



V° PE (p',p;E) 



2m 



J dzg(z, p')f(p', p, z)g(z, p)g(E - z, -p') 



J dz g(z, p)g(E -u k -z, -p') J dz g(z, -p')/(-p', -p, z)g(z, -p)g(E - z, p) 

(50) 

The fully off-shell OPE potential V^ PE (p', p; E) is related to the OPE Green function 
in the usual way: 

V° PE (p', p; E) = G^(E, p\ - p ')V^ p V, P; E)Gnn(E, P, -p) (51) 



where Gnn is the dressed two- nucleon propagator given as in Eq. (|24[ ) . For the numer- 
ical comparison of Eqs. (HHf) and (|50) we follow our previous work |L3| and use the Ml 



ttN interaction of Ref. JTTJ. We shall, in each case, calculate the half-off-shell poten- 



tial V^ PE (p , p; E), defined by Eq. (|5T| ) in the limiting case p' — > p , where p is the 
on-shell momentum, i.e. E = p^/m + 2m. Here we shall consider the simplified case 
where no dressing is included for the vertices of the exchanged pion; in this case, both 



the two-energy vertex f(p' } p, E', E) of Eq. (|48|) and the single-energy vertex f(p', p, E) 



of Eq. (|5l]|), reduce down to the energy-independent bare vertex /o(p')P)- Note that the 
restriction to bare vertices is limited to the pion exchanged between the two nucleons, 
so that the dressed nucleon propagators g are totally unaffected by this simplification. 
The half-off-shell potential V^ PE (po, p; E) depends on three variables, the energy E, the 
magnitude p of the off-shell momentum, and the cosine of the angle between po and p. 



Figure 4: Comparison of the half-off-shell dressed NN OPE potential, V^ PE (p , p; E) 
calculated with full dressing (solid curves), and within the convolution irNN model (long- 
dashed curves). The short-dashed curve is for the case where no dressing is included. 

For the numerical comparison, we have examined the potential as a function of energy 
E for a large range of values of p and x = po ■ p. A typical result is shown in Fig. |4] 
where, in this case, we have set p = .lpo and x = —0.8. In this figure, the solid curve 
gives the OPE potential with full dressing, while the long-dashed curve is the result when 
connected three-body forces are neglected. Also shown is the standard OPE potential 
where no dressing at all is included (short-dashes). The first observation of note is the 
significant effect that dressing has on the standard OPE potential. This raises questions 
about the role such dressing may play in standard descriptions of the NN force. The 
second observation forms the main result of this section - it is the essential identity of the 
solid and long-dashed curves, indicating that the contribution of connected three-body 
forces to the OPE potential is negligible. Although strictly applying to the OPE poten- 
tial with bare vertices, this is a very encouraging result that gives us hope that connected 
three-body forces may be small in general. If this is borne out in further studies, this 
would mean that, for the ttNN system, the convolution ttNN equations provide a way of 
effectively summing all possible diagrams of time-ordered perturbation theory. 

DRESSING OF CONNECTED DIAGRAMS - METHOD I 

In Ref. [|nj we showed how a convolution integral can sum all possible relative time 
orderings of a disconnected diagram of time-ordered perturbation theory. This was then 
used to demonstrate how two disconnected nucleons can be dressed simultaneously. Here 
we show that the convolution idea can also be easily extended to connected diagrams. 
The final result shows that the sum of all topologically similar diagrams, connected or 
disconnected, differing only in the relative time ordering of their vertices, can be expressed 
through a convolution formula that effectively integrates out initial and final relative 
energies from the topologically equivalent Feynman diagram. One application of this 
generalized convolution formula is to the problem of including all possible dressings in 
connected diagrams. 

The starting point of the following discussion is relativistic quantum field theory. Since 
Feynman diagrams contain, in some sense, all time orderings, it is not surprising that our 
goal will be to show how to extract the Green function of time-ordered perturbation theory 
from the Green function of relativistic quantum field theory. Although the procedure we'd 
like to follow is general, for presentation purposes we specifically consider the NN — > NN 
process where the interaction is described by a Hamiltonian H = Hq+Hj involving meson 



and baryon fields. The explicit form of H need not be specified. 

The free fermion field at time t = is denoted by ^>(x). At time t, we define 

V>(x,£) = e iHot ip(x.)e- iHot 
tf(x,t) = e iHt ip{x)e- iHt 



(52) 
(53) 



being the interaction picture and Heisenberg fields, respectively. The free fermion field 



(54) 



tjj(x) = ip("x., i) can then be written in terms of its Fourier decomposition as |L8| 

We consider the NN^-NN process described by the coordinate space Green function 

ig(x,y;x',y') = ((0|r*(a:)*(i/)*(a/)*(j/)|0», (55) 

which, because each \l/ is a four- component spinor, can be considered as a 16 x 16 matrix. 
Here |0)) is the dressed vacuum, its relation to the bare vacuum |0) being given by 



|0» 



1 [l+^ff, 



((0|0) L ~ ' 0+-H 
The momentum space Green function is defined by 

i(x-p+yq-x'-p'-y'-q') 



G(p,q;p',q 



Q(x,y; x' \y')d A x d A y d A x' d A y' 



(56) 



(57) 



which can also be expressed without the momentum conserving 5-function by defining the 
Green function G(p, q; p', q')\ 

G(p,q\p',q') = (2tt) 4 5V + q' - p - q)G(p,q;p',q'). (58) 
We follow Logunov and Tavkhelidze [T9] and consider the two-time Green function 



C?(x, y, t; x, y', t) = Q(x, y;xy 



xq = yo = t ■ 



(59) 



The two-time Green function in momentum space is defined by 
g(p,q,E;p',q!,E') = J e H-^P-y^W-p'+yW+tE-fE') 

(?(x, y, t; x', y', t')d 3 x d 3 y d 3 x' d 3 y' dt dt' . 



(60) 



A little algebra shows that 



C?(p, q> E; p'> q'> 



l 



(27r)^ J —oo J— oo 



dull duj' G(p, E — u, q,cj;p', E' — to' , q' ',&/)■ (61) 



This result was used by Logunov and Tavkhelidze as the starting point for their study of 
the two-time Green function. It also constitutes the solution to our problem as we now 
proceed to show. 

It is useful to express Q without momentum conserving (^-functions, thus we define the 
Green function G by 



0(p, q, E] p', q', E') = (2n) 4 5(p + q - p' - c[)5{E - E')G(E, p, q; p', q'). (62) 



Removing the momentum conserving ^-functions from Eq. (|6l|) , gives 

1 



G(£,p,q;p',q' 



du 



(271")*^ J—oo J — oo 



du' G(p, — u, q, cu; p', E — oj', q', a;'). (63) 



The rhs of Eq. (|63|) is a double convolution over initial and final relative energies 
of Feynman graphs. The goal, therefore, is to show how the lhs, G(E, p, q; p', q'), is 
related to the Green function of time-ordered perturbation theory. For this we go back 
to coordinate space and examine the two-time Green function in detail. By the definition 
of the time-ordered product in Eq. (|55|), 

^(x,y,t;x',y',t') = 9(t - *')«0|*(x, *)*(y, t)*(x', t')\0)) 

+ 9(t> - t) «0|*(x', t'My', 0*(x, f)*(y, t) |0». (64) 

Using Eq. (|53"D in this equation, and taking the Fourier transform as in Eq. (|60"D gives the 
result 



(27r) 3 5(p + q - p' - q!)G(E, p, q; p', q' 



((0|V(p)^(q 



1 



where we have used Eq. (|62|) , the fact that 

9(t)e iujt dt 



lu + ie 

and where we have introduced the momentum space fields 



E+-H 

PW)^^^(pMq)|0>> (65) 



(66) 



V>(p) = / rfx 3 e- ix - p ^(x). 



(67) 



Using Eq. (|>|) at t = in Eq. (|67D , one obtains that 



^(p) = E 



(2vr) 3 m 



[6(p, s)u(p, s) + d f (-p, s)v(-p, s)]. 



(6? 



Eq. ( |65"D is basically the result that we seek. It expresses G in terms of two terms corre- 
sponding to the retarded and advanced parts of the Green function of Eq. flB3|); we shall 
correspondingly name the first and second terms on the rhs of Eq. (|65|) as retarded and 
advanced, respectively. 

One may compare Eq. (|65"D with the Green function of time-ordered perturbation 
theory, which for the process in question is given by 



(2vr) 3 (5(p + q - p' - q')G a(3 ^ p/ (E, p, q; p', q') = (p, a, q, (3\ 



E+-H 



pV,q',/?'> (69) 



where the spin components a, (3, a', j3' are shown explicitly. Taking into account Eq. 
we see that the retarded part of G in Eq. ( |65|) is very similar to the Green function of Eq. 
(|69|). This similarity suggests a simple transformation of Eq. (R3) defined by the equation 



G 



a/3, a' /3' 



(27r) 3 m 



EpE q u(p, a)w(q, (3) 



G 



v Ep,E q , 
{1-K) 6 m v 



(70) 



where here G = G(p, q, p', q') represents any appropriate Green function. Note that Eq. 
(|T0| ) transforms a 16 x 16 matrix in spinor space to a 4 x 4 matrix in spin space. Under 
this transformation, Eq. ( |65f) becomes 

(2vr) 3 5(p + q - p' - q')G a/W (£, p, q; p', q') = 

«0|6(p, a)6(q, fl^J-^V, a W, /3')|0» 

- (<O|6t(,/ I a06 t (q / ,^^ ¥ 6(p,a)6(q I /3)|O>>. (71) 



If |0)) = |0), the advanced term of Eq. (|7TJ) disappears, and the retarded term becomes 
identical with the Green function of time ordered perturbation theory, Eq. (|69|) . Thus 
the only difference between the transformed two-time Green function of Eq. fl7Ip and 
the one of time-ordered perturbation theory, is in the type of vacua used to take the 
matrix element: the transformed version of G uses the dressed vacuum while standard 
time-ordered perturbation theory, as in Eq. (|69|), uses the bare vacuum. 

In general, |0)) ^ |0), however, in the special case where the dressed and bare vacua are 
equal, Eqs. ( j69| ) and ([H]) become identical. In this case, Eq. (|63]), after the transformation 
of Eq. (f70D, provides us with a formula that expresses the full Green function of time- 
ordered perturbation theory in terms of a double convolution integral of the full Feynman 
Green function. This result can be easily generalized to hold for perturbation diagrams 
of a given order in the interaction, or even for perturbation diagrams of a given topology 
(the argument is similar to the one used in Ref. fl3fl). It is in this latter form that the 
method can give particularly useful formulas, like that of Eq. (|48|) , where all possible 
time-orderings of a perturbation graph are summed by performing convolution integrals 
of one Feynman diagram. 



A similar result to ours has lately been obtained by Phillips and Afnan J2C] using 
a much more involved argument. In their case, however, they neglected antinucleons 
altogether in order to obtain the connection of Eq. fl6"3|) to time-ordered perturbation 
theory.0 

We emphasize that all that is needed in our derivation is the condition |0)) = |0). 
For this to be true, it is sufficient to drop terms from the interaction Hamiltonian that 
connect vacuum to vacuum states, otherwise antinucleons can be retained. For example, 
in the simple case of interaction Hi ~ 4>4>ip, expanding each field in terms of creation and 
annihilation operators gives eight terms; however, of these eight, only terms b^a^S and dab 
contribute to the dressing of the vacuum (in lowest order via the process — > NNn — > 0), 
and just these terms can be dropped while still retaining other terms involving d and d) . 



DRESSING OF CONNECTED DIAGRAMS - METHOD II 

The above method using the two-time Green function can be used to derive Eq. 
for the dressed AW OPE potential, but only in the case where vacuum dressing has 
been neglected. It turns out, however, that Eq. (|48|) is more general, holding for any 
Hamiltonian H, including ones involving vacuum dressing. 

Here we derive Eq. (IOT) using a totally different method that is close in spirit to the 



one used to derive the convolution formula for disconnected diagrams | 13| . This method 



2 This paper claims, incorrectly, to have shown that the convolution formula, derived by us in Ref. 
[fl3|| , is a special case of their two-time approach where antinucleons are neglected. If anything, just the 
opposite is the case: our convolution formula was derived using a method, akin to Method II, which does 
not need for its validity any approximations whatsoever. We also do not share these authors' criticism 
of three-dimensional approaches, even though we have ourselves already presented the first consistent 
four-dimensional approach Q. 



involves the temporary introduction of different types of pions and nucleons into a time- 
ordered perturbation theory description, and provides a way to derive Eq. fl4"B|) without 
making any approximations. In this sense, this method is more general than the one 
based on the two-time Green function. 

The central idea is to interpret the dressed NN OPE potential as consisting of two 
distinguishable nucleons Ni and N 2 , each being dressed by its own pion, tti and ir 2 re- 
spectively, and in addition, exchanging a third type of pion 7r; of the three types of pion, 
only 7r can interact with both nucleons. We accordingly define Hi = H (l) + Hj(l) to 
be the Hamiltonian describing the ttiNi system, and similarly H 2 = Hq(2) + Hi (2) is 
the Hamiltonian describing the ir 2 N 2 system. The Hamiltonian describing pion 7r and its 
interactions with the nucleons can be similarly written as H n = Hq + Hf(l) + Hf(2). 

Thus, given any Hamiltonian H = H Q + Hj for which we would like to calculate the 
NN OPE potential, we can proceed by replacing H with the sum H' = H\ + H 2 + H 71 , 
where each individual Hamiltonian Hi, H 2 , and H n has free and interaction parts that 
are of the same form as H and Hi, respectively, but each with its own individual fields 
replacing the corresponding ones of H. 

Although it may be possible to provide a general formulation, it is more convenient 
to illustrate the procedure to follow by taking the usual model where the interaction is 
given by a three-point ttNN vertex. 

As we shall not explicitly need to use any details of the model for Hi and H2, we 
specify the model in terms of the Hamiltonian involving the pion 71: 



H£ = J dkw fc 4(kK(k) (72) 
HJ = [ dkal(k)J N (k) + H.c. (73) 
J N (k) = [dpdp'6(p + k-p')-^F (p,p')a ] N (p)a N (p') (74) 



where iV can be either Ni or A^2, in which case HJ needs also to be labelled accordingly. 
Note the relations 

Mk),#J] = co k a n (k) (75) 
[a w (k),HJ] = J N (k). (76) 

Consider now the perturbation expansion of the full NN —>■ NN Green function with 
respect to the interactions Hf(l) and H]{2): 

(Pi,P^i^r^7lPi,p 2 ) = (p^ V-^-^-ff^ 1 ) 

E+-Hi-H 2 - HZ HJ{2) E+-Hi-H 2 - H S ,Pl ' P2) + ' ' ' (77) 

where the term of order Hf(l)Hf(2) has been singled out of the complete perturbation 
series, as it is just this term that coincides with the exact OPE potential specified by 
the original Hamiltonian H. Note how the introduction of the three Hamiltonians Hi, 
H 2 , and H*, enables us to treat meson exchange perturbatively, while nucleon dressing is 
treated non-perturbatively. 

Now consider only this OPE term. Replacing Hf(l) and HJ(2) by the integral of Eq. 



(|73|), and then using Eq. ([75|) , we obtain that 



G ?r = / rfk (P ;, rt | E+ _ Hi _ H2 jU ( k) E+ _ Hi _ H2 _ Uk j^m 

1 



where we have used that ifg acting on two-nucleon states gives zero, and where one k'- 
integral has been eliminated using [a w (k'), aj.(k)] = 5(k' — k). The essential step comes 
at this stage when we recognize that [Hi, H 2 ] = [Hi, J^ 2 ] = [H 2 , JjvJ = which enables 
us to write Eq. ( |78"D in terms of two contour integrals: 

1 . + , 1 , . 1 



Hi JV1V z + — Hi z' — z — ujk 

where the matrix element factors into two, one factor for each nucleon. Because of this 
factorization, and because all Hamiltonians have the same form as H, we may now drop 
the nucleon labels in Eq. (|79|). 

The matrix elements in Eq. (|79|) define the two-energy vertices 

6(p' + k - p)/(p' ,p, z' , z) =g~ 1 (z',p) z ,+ \ h J n{^) z+ \ H \v) g~ l {ziV>)- (80) 

Substituting this definition into Eq. (|79|) results in the expression of Eq. (p8|) . 
Using Eqs. ( [751) an d (|7B|), it is straightforward to show that 

<p'i?^V»(k)^ip> 

= (p ' k| ?~ff |p> + {z ~ z ' ~ " Jk){p ' l 7^H a ' {k, 7^ L H lp) ■ (81) 

Recognizing that the first term on the rhs is the usual one-energy vertex function, Eq. 
(|8Tp shows that the one- and two-energy vertex functions coincide for on-mass-shell pions. 



SUMMARY 

Attempts to formulate few-body equations for the ttNN system have a long history. 
Yet all attempts from the very first bound state model of the 1960 's to the sophisticated 
NN — ttNN models of the 1990's had one basic feature in common, each attempted to 
simplify the field theory in question by truncating to some maximum number of pions. 
Unfortunately, this truncation results in an inconsistent treatment of nucleon dressing with 
serious consequences for practical calculations. The NN — ttNN equations, for example, 
suffer from a renormalization problem which effectively reduces the strength of the input 
ttNN vertices, which in turn contributes to the underestimation of cross sections. 

In order to resolve these renormalization problems, we have presented a completely 
different approach where the Hilbert space is not truncated to some maximum number 
of pions. Instead, the guiding principle has been to work only with fully dressed vertices 
and propagators. What this means in practice, is that all possible disconnected diagrams 
of quantum field theory have been retained in our model. To achieve this, we have used 
convolution integrals to sum over all relative time orderings of disconnected graphs. As 
convolution integrals have often been used in nuclear physics, especially in the four-nucleon 
problem ^TJ, it may appear surprising why the idea has not been applied sooner to the 



7riV7Vproblem. This, however, can be understood when it is realized that the convolutions 
that we have used are applied to field theory Hamiltonians, and result in convolution 
expressions for energy- dependent" potentials" . This is quite different from the more usual 
convolutions of nuclear physics which apply only for energy-independent potentials. 

The convolution irNN equations are unitary and do not suffer from renormalization 
problems. These are the two crucial attributes that are necessary for a theoretically 



consistent description. Indeed, we should emphasize that the downfall of the NN— ttNN 
equations is not because the two-nucleon propagator is underdressed, rather, it is simply 
because the equations themselves do not have both the properties of unitarity and correct 
renormalization. 

The only approximation made in deriving the convolution ttNN equations is the neglect 
of connected three-body forces. This is, of course, also an assumption in essentially every 
model in nuclear physics. Drawing on the bulk experience of the field, such three-body 
forces are very likely to be small. However, to obtain a more quantitative assessment, 
explicit calculations of connected three-body force contributions are needed. With this 
goal in mind, we have extended the convolution idea also to connected diagrams. 

Two methods have been presented to derive convolution integrals that sum all the 
time-orderings of a connected graph. The first method, involving the two-time Green 
function of relativistic field theory, is useful to obtain the form of the convolution integral 
in a quick and straightforward manner. However, this method works only in the case when 
dressed and bare vaccua are equal. The second method, is an extension of the one used by 



us to derive the convolution formula for disconnected diagrams ||13||, and involves replacing 



the true Hamiltonian H by a sum of Hamiltonians, each having the same form as H, but 
applying only to singled out particles. This method works without any assumptions and 
is therefore more general than the one involving the two-time Green function. 

With the convolution formula for connected diagrams derived, we have calculated the 
fully dressed NN OPE potential, since part of this potential includes connected three- 
body forces. For this dressed OPE potential, we find that the corresponding connected 
three-body force is negligible. Thus we can conclude, that at this stage, the convolution 
ttNN model is a candidate for effectively summing the whole of time-ordered perturbation 
theory for the irNN system. 
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